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,^ ' A semiclassical nonlinear collisional drift wave model for dense magnetized plasmas is developed 

.^- 

and solved numerically. The effects of fluid electron density fluctuations associated with quantum 

►-j>. , statistical pressure and quantum Bohm force are included, and their influences on the collisional 

00 

f~^ \ drift wave instability and the resulting fully developed nanoscale drift wave turbulence are dis- 

cn ■ 

cussed. It is found that the quantum effects increase the growth rate of the collisional drift wave 
instability, and introduce a finite de Broglie length screening on the drift wave turbulent density 
perturbations. The relevance to nanoscale turbulence in nonuniform dense magnetoplasmas is 
discussed. 
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Recently, there has been considerable interest \Mi^ in the investigation of collective dy- 
namical processes in dense quantum plasmas, which are ubiquitous in astrophysics (e.g. 
white dwarf stars), in planetary systems (e.g. Jupiter and giant planets around extraterres- 
trial stars), in plasma-assisted nanotechnology (e.g. quantum diodes, quantum free-electron 
lasers, metallic thin films and nanostructures, nanowires, nanoplasmonics, etc.), and in 
highly compressed plasmas for inertial confinement fusion by extremely powerful laser and 
charged particle beams. 

In quantum plasmas, the dynamics of degenerate electrons is governed by quantum hy- 
drodynamic (QHD) equations [l| in which the electron momentum equation has novel forces 
associated with the electron pressure involving the Fermi-Dirac distribution function, and 
the quantum Bohm force arising from the overlapping of electron wavefunctions due to 
Heisenberg's uncertainty principle. Furthermore, in a dense plasma non-degenerate ions are 
in a strongly coupled state since the ratio between the Coulomb ion interaction energy and 
the ion kinetic energy is much larger than unity. Thus, kinematic ion viscosity and vis- 
coelastic ion relaxation for ion correlations have to be included in the description of the ion 

nn 

motion. It then turns out that collective effects pj-|^ in quantum plasmas arise at nanoscales 
due to relatively high plasma number density and the low electron and ion temperatures 
when compared with the classical plasma. We thus have new quantum regimes, the physics 
of which is quite interesting and appealing for practical applicatioris, as mentioned above. 

The excitation of the electrostatic (ES) drift waves (DWs) p, |6| is a universal feature 
of the classical nonuniform magnetoplasma. Classical ES DW turbulence has been mostly 
investigated by using the two- fluid theory and simulations [3, 18[, where non-degenerate elec- 
tron and ion fluids in the presence of the low-frequency (in comparison with the ion gyrofre- 
quency) drift waves are dynamically coupled via self-consistent ES fields and finite density 
fluctuations. It is widely thought that DWs are responsible for the cross-field transport of 
plasma particles in magnetic fusion devices J7|. The two- fluid quantum magnetohydrody- 
namic (Q-MHD) equations for a uniform quantum magnetoplasma, by including electron-1/2 
spin effect, have been developed by Brodin [9|, by assuming that non-degenerate ions are 
uncorrelated. Furthermore, the linear properties of the ES DWS in a nonuniform quantum 



magnetoplasma have been studied by several authors [id. Ill[|. 

In this Letter, we develop nonlinear equations for the low-frequency ES DWs in a nonuni- 
form quantum magnetoplasma that is collisional. The governing mode-coupling equations 



are then numerically solved to depict the novel features of fully developed ES DW turbu- 
lence at nanoscales in a driven (due to the combined action of the density gradient and 
electron-ion collisions), nonuniform Fermi magnetoplasma. 

Let us consider a nonuniform quantum magnetoplasma in an external magnetic field B. A 
background plasma density gradient Vno(x) is perpendicular to the magnetic field direction 
b = B/S along z in a local Cartesian coordinate system. 

The governing nonlinear fluid equations for electron with density ?T,e(x, t) and velocity 
v(x, t), and ions with density nj(x, t) and velocity u(x, t) are the continuity equations 

(9i + V ■ V)ne + neV ■ V = 0, (1) 

{dt + n- V)ni + UiV ■ u = 0, (2) 

and the momentum equations. The inertialess momentum equation for degenerate electrons 
is 

Fe Vpe + FB-meI/ei(v-u) =0, (3) 

Tig 

where Fg = e (V0 — v x B/c) is the Lorentz force on electrons with charge e in a dynamical 
electrostatic potential 0(x, t) and static magnetic field B. 

The pressure term for the non-relativistic degenerate electrons is expressed as pe = [Te + 
(5/3)Tp]ne, where the thermal electron Temperature Tg is in the following assumed to be 
much smaller than the Fermi electron temperature Tp = (/i^/2r?7,e)(37r^?T,o)^^'^. The Fermi 
energy for electrons with mass rrif, is rngVp where Vp = \/Tpjmf^ is the Fermi velocity at 
equilibrium density no. 

The effect of quantum hydrodynamical diffraction on electrons is included by the Bohm 
force F^ = {7i j^mf^SJ (V^ In y^) . Electrons are coupled to ions via the collision frequency 
z/ei representing the momentum exchange, with mjnjZ/je = ■me^e^'ei- 

The momentum equation for inertial ions is 

ramiidt + u ■ V)u + n^i + Vpj + njmji/je(u-v) = R (4) 

where the viscous term for strongly coupled ions is determined by 

(1 + r^9i)R = r/V- Vu + (e + r7/3)V(V-u) (5) 

with viscoelastic relaxation time r^, and ri and ^ are the longitudinal and bulk viscosities 
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13|, respectively. 



The Lorentz force on ions with charge state Z is Fi = Ze (V0 — u x B/c). 

The ion pressure is pi = •yirii^fiiTi + T,,,), where 7^ is the ion adiabatic index. The excess 
chemical potential /ij ~ 1 + 0.33U(Ti) + (Ti/9)dU(Ti)/dTi, includes U(Ti) as a measure of 
the excess internal ion energy, where Fj = Zf/aiTi, ai = (S/ATrriiY^^ is the Wigner-Seitz 
ion radius, and Tj the ion tem per ature. T^ = {Nn/3){Zfe'^/ai){l + K)exp(— k) accounts for 
the strong ion coupling effect [ij]. Nn is determined by the ion structures and corresponds 
to nearest neighbors ion interactions if the ions are in a crystalline state, k = ai/XTF, and 
Xtf = vp/ijJpe the Thomas-Fermi screening radius, with ojpe = {A-nriee^ /meY^"^ being the 
electron plasma frequency. 

In the following the drift approximation is applied, assuming that all dynamical frequen- 
cies u <^ Uce are low compared to the electron gyro frequency u^e = eB/rrieC. The velocities 
are expressed as v = v^ + bt>||. The perpendicular component of eq. [3] gives the electron 
fluid drift velocities v_|_ = v^; + vp^ + Vg, where v^; = {c/B"^) B x V0 is the classical 
E X B drift velocity, \f* = {cTp / ennlB^) B x Vn^ is the electron Fermi drift velocity, 
and Vg = (ch /2eB'^men) B x V_l (V^ In y/n) is the quantum diffraction drift velocity. Here 
we have assumed that Uei ^ <-^ce- Quasi-neutrality requires he ^ hi = n. The electron 
continuity equation in drift approximation is written as 

{dt + VE-V)n + nV-VE-V ■ {nvp*) - V ■ (nvg) - V||H^) = 0. (6) 

In addition to the dominant convecting E x B velocity, the drift velocities acquire a finite 
divergence in an inhomogeneous magnetic field, which enters like a compressibility into the 
continuity equations by abbreviating the curvature operator as 

Kif) = -cV X (B/52) . V/ = -cV X [(B X V/)/i?2] . (7) 

We obtain the divergence of the drift velocities as V ■ v^ = —k{(J)), V ■ (nvi?,,,) = (TF/e)K{n), 
and V ■ (nvq) = —(fi K/Amee)V'^n. The density n = no{x) + h{x,y,t) can be now split into 
the background and fluctuating components, viz. uq^x) and h{x,y,t), respectively, where 
h{x,y.t) <ti no{x). 

The background gradient advection term can then be expressed as v^; • Vno = 
— {c/LnB)dy(j), with L~^ = \dx\nno\. In the following we drop the tilde on the fluctuat- 
ing density component for clarity. Equation (6) then becomes 

Dtn = ^dy4> - -ViiJii - ^K{n) + K{<t>) - -^K (V^n) , (8) 

BL„ e e Arripe ^ 



where Dt = 9t + vg ■ V. The parallel electron current density j\\ is obtained from the 
momentum equation as 

i|| = ^^V||0-^^V||(l-AjV> (9) 



where the electron de Broglie length at the Fermi temperature is Xq = %/ y/Am^Tp. The ions 
are assumed two-dimensional with f|| ^ u\\ = 0. 

For non-degenerate ions, we assume that Tmdtu± ^ u^, where the perpendicular com- 
ponent of the ion fluid velocity in the drift approximation {u ^ Ud = ZieB/rriic) is 
u± = v^; + Up. The E X B drift velocity for ions u^; = v^; is identical to those for 
electrons, and the ion polarization drift velocity is Up ^ —{c/ BuJci)DtV ^(j). Here the 
ion gyrofrequency Ud is assumed to be much larger than /Xje, and we have assumed that 
l-D^I ^ [{i + 4?7/3)/nomj]V5_- Assuming Ti^ T^ <^ Tp, we have also excluded the contribu- 
tion of the ion diamagnetic drift in u^. The ion continuity equation in drift approximation 
then yields 

Dtn = ^dyd^ + /.(0) - 1^ AVl0, (10) 

Subtracting (IHl) from ( TTOl) we obtain the modified ion vorticity equation 

= AVi-f ^Viuii - ^-^^K{n) - ^/t V^n . 11 

ip riQec HqCC AenorrieC 

Eqs. dH]) and (TTTj) constitute coupled nonlinear dynamical equations for low-frequency 
perturbations of the density n and electrostatic potential (p in an inhomogeneous magne- 
tised semi-classical plasma. The standard drift normalization is applied, with the Fermi 
temperature Tp as reference. Spatial scales x/Lj_ — )■ x are in units of a perpendicular refer- 
ence length L±, and times scale as tCs/L± — > t with Cs = ^/Tp/mi. Fluctuating components 
are normalized as S^'^ecp/Tp — )■ and 5Q^n/no — )■ n, where ^o = p/L± with Fermi drift scale 



p = y/Tpmi/{eB). The gradient length scale enters as gn = L±/Ln. Two possibilities for 
length scale normalization are appropriate, setting L± either to p or to L„. 

The parallel derivative acting on the current is not explicitly evaluated here, but rather 
determined by a fixed parallel wavenumber k\\ through a dissipative coupling parameter 
d = {LOce/i^eiSo)kfi- Accordingly, we can rewrite eqs. ([8]) and ( TTT]) as a set of normalized 
quasi-two-dimensional equations 

dtn + [0, n] = d{(t) - An) - K{K*n) (12) 

dtn + [0, n] = d{(f) — An) — gndycf) + K{(f) — A*n) (13) 



where fl = V^0 and [a, b] = z x Va ■ Vb represent the ion vorticity and Poisson's bracket, 
respectively. The quantum diffraction effect enters through A = 1-/3^ V^ and A* = 1+/3^V^, 
with /3 = Xq/p. It acts as a finite de Broghe length screening effect on density ffuctuations, 
similarly to the established FLR gyro screening and Debye screening effects. Eqs. flT^ and 
[T5]l are the semi-classical dense plasma generalization of the Hasegawa-Wakatani equations 
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First, we derive the linear dispersion relation from Eqs. f lT2|) and (1131) for k = 0. For 
this purpose, we neglect nonlinear terms and suppose that and n are proportional to 
exp(— zwt + zk-r), where u and k are the normalized frequency and wave vector, respectively. 
We have 

k'^uj'^ + id{l + Ak'^)Lo - idAgnky = (14) 

where A = (1 + /3^A;^). In the adiabatic limit d ^ 1 we obtain w* = Agnky/{1 + Ak'^) as the 
electron diamagnetic drift frequency including quantum corrections. Inserting u = ujr + ^70 
into the dispersion relation and solving for the imaginary component in the limit 70 ^ w* 
we obtain the linear growth rate for weakly nonadiabatic quantum drift waves: 

^° - rf(l + AP)3 = ^ ^^^> 



with Wo = d{l + Ak^)/{2k^). The exact solution 

V2 



7 = Wo 



1± 
2 



y^l + Vl + 16(7oM)' 



(16) 



for the imaginary part of the frequency from dispersion relation flT^ is shown in Fig. [T] for 
kx = 0, On = ^, d = 0.1 (black) or d = 2 (red), and various values of /3. The maximum 
growth is found around pky ~ 1. In the more strongly coUisional case {d = 0.1), the 
maximum growth rate is increased by around one third for /3 = 1 compared to /3 = and 
is shifted towards slightly smaller wavelengths. For lower coUisionality {d = 2), the growth 
rate remains nearly constant when /3 is raised, but the maximum shifts to slightly lower ky, 
while the de Broglie screening significantly reduces the growth rates for pky > 1. For d > 2 
the growth rate 7 ~ 70 is well approximated by the weakly adiabatic solution. 

Next, we numerically solve the nonlinear quantum drift wave Eqs. (1121) and (I13p for 
a uniform magnetic field {k = 0). For time stepping, an explicit third order Karniadakis 



scheme [l6| is applied, and the Poisson bracket [a,b] = {dxa){dyb) — {dya){dxb) is evaluated 
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FIG. 1: The linear growth rate j{pky) for d = 0.1 (black) and d = 2 (red), with the de Brogl 
factors (3 = (hold curves), (3 = 0.2 (thin curves) and (3 = 1.0 (dashed curves). 



with the energy and enstrophy conserving Arakawa method p/7[. The numerical method 
is equivalent to the one introduced in Ref. [18[. A hyper- viscous operator, e.g. associated 
with strong ion coupling effects, z/^V^ with z/^ = —10^^ is added for numerical stability 
into the right-hand side of both Eqs. [12] and [131 acting on Q and n, respectively. The 
Poisson equation is solved spectrally. The equations are for the following computations 
discretized on a doubly periodic 1024 x 1024 grid with box-dimension L = 64p. Spatial 
scales are in units of p. Here we choose L± = Ln and thereby gain fluctuations and n 
in the order of unity (for 6o ^ I and e^t/Tp ^ 1, u/uq <^ 1)) and accordingly set Qn = 1. 
The nominal parallel coupling parameter is set to d = 0.5. The computations are initialized 
with a Gaussian density perturbation which transiently develops into a drift wave during a 
quasilinear instability phase and finally saturates nonlinearly into a fully developed turbulent 
state. 

In Fig. [2]the wavenumber spectra of the electrostatic potential fluctuations |0(A;)| (black) 
and density fluctuations \n(k)\ (red)in the fully developed turbulent state are shown for 
/3 = (bold) and 0.5 (thin). The spectral properties of |0(fc)| remain similar with larger /3, 
while its fluctuation level is increased for finite /3 for all wavelengths as a consequence of the 
larger growth rates. The density |^(/i:)| on the other hand shows a strongly reduced level 
for all but the largest scales. For a larger dissipative coupling parameter d (corresponding 
to a lower electron-ion collision frequency) the electrostatic potential fluctuations are more 




FIG. 2: Spectra \(j){k)\'^ (black) and \n{k)\'^ (red) of a fully developed quantum drift wave turbulence 
for (3 = (bold) and 0.5 (thin). 



strongly coupled to the density perturbations and experience similar screening. 

Our simulations show for the first time the effect of de Broglie screening on dense plasma 
turbulence. This novel quantum diffractive screening mechanism also affects the turbulent 
fluctuation energies. The total energy E = (1/2) J dV (ri^ + (V0)^) is only slightly changed 
(by around ±10%, depending on the plasma parameters) for (3 < 0.5 compared to /3 = 0. 
Its components, the potential energy En = (1/2) J dVn^ and the kinetic energy Ek = 
(1/2) J dV('V(f)y), are however strongly and reciprocally influenced. Figure [3] shows En{t) 
(top) and Ek(t) (bottom) for /3 = (bold) and 0.5 (thin). While the density perturbations 
decrease with increasing /3, the velocity fluctuations increase. The quantum effect on the 
drift wave turbulence is visualised in 2-D plots of the turbulent fluctuations in Fig.lHfor /3 = 
(left column) and 0.5 (right column): the de Broglie screening affects density perturbations 
(top row), but conserves the fine structure in the ion vorticity (bottom row). 

In conclusion, we have derived and solved nonlinear DW equations for nonuniform Fermi 
magnetoplasmas that are collisional. The present nonlinear equations include the new 
physics of quantum forces involving the pressure of degenerate electron fluids and over- 
lapping electron wavefunctions over nanoscales. The quantum regime is characterized by 
P = Xq/p, which is the ratio between the de Broglie length and the drift scale at the 
Fermi electron temperature. The quantum effect acts as a de Broglie screening on density 
perturbations for small wavelengths, and increases the growth rate and turbulent velocity 
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FIG. 3: Time traces of the potential energy r? (top) and the kinetic energy {Vcl))'^ (bottom) for 
/? = (bold) and 0.5 (thin) 



fluctuations across the whole spectrum. The DW vortex structures occur at nanoscales, 
which are of the the order of several Fermi sound gyroradius p. The novel de Broglie effect 
on turbulent drift vortices is of a similar quality as finite Larmor radius (FLR) corrections 
in gyrofluid or gyrokinetic simulations of warm ion (ITG) drift wave turbulence [l9J, and as 
Debye shielding effects on electron gyro scale (ETG) turbulence [20i] . 

In closing, the results of the present investigation were primarily motivated by theoretical 
interest, but may be useful in understanding turbulence and coherent structures that arise 
in nonuniform dense magnetoplasmas, such as those in the atmospheres of white dwarfs 
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241 ] ■ Fully developed DW turbulence can regulate the cross-field plasma particle transport 
over nanoscales. 

For example, a mixing length estimate of DW turbulent convective particle diffusiv- 
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uj^:p{TpY in comparison with collisional diffusive cross-field transport D, 
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UeiPiiTiY in a strongly magnetized plasma depends essentially on the the (local or radially 
global) density gradient length, which is subject to a large uncertainty. Assuming a white 
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FIG. 4: 2-D plots of the density n{x, y) (top) and the ion vorticity Q(x, y) (bottom) for /3 = (left) 
and (3 = 0.5 (right). Only a quarter section ('512 x 512j of the computational domain is shown. 

dwarf surface density n ~ 10^°c?7i~^, the coUisionality is in the order of Uei ~ lO^^s^^ (in- 
dependent of density in a degenerate plasma). The Fermi sound speed is in the order of 
Cs ~ lO'^cm/s (smaller than 0.1 c, justifying the present non-relativistic approximation), and 
the radially global density gradient length is in the order of L~^ ~ lO^cm"^, giving w* = 
c,/L„ ~ lO^s'S so that i/^iM ~ 10^. On the other hand, piTpf / pi{T,f = Tp/Ti ~ 10^ 
is of the same ratio (for an ion temperature Tj in the order of lO^i^), so that Dt^rb ~ Dcoii 
are in the same range. Passive cross-field advection of heavier trace ions by small-scale 
DW turbulence could contribute to counter sedimentation in cool white dwarfs, even in the 
absence of thermal convection or radiative acceleration. 
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